Phase Transitions in a Kinetic Flocking Model of Cucker-Smale Type

We consider a collective behavior model in which individuals try to imitate each others' velocity and have a preferred speed. We show that a phase change phenomenon takes place as diffusion decreases, bringing the system from a “disordered” to an “ordered” state. This effect is related to recently noticed phenomena for the diffusive Vicsek model. We also carry out numerical simulations of the system and give further details on the phase transition

Propagation of chaos for rank-based interacting diffusions and long time behaviour of a scalar quasilinear parabolic equation

We study a quasilinear parabolic Cauchy problem with a cumulative distribution function on the real line as an initial condition. We call 'probabilistic solution' a weak solution which remains a cumulative distribution function at all times. We prove the uniqueness of such a solution and we deduce the existence from a propagation of chaos result on a system of scalar diffusion processes, the interactions of which only depend on their ranking. We then investigate the long time behaviour of the solution. Using a probabilistic argument and under weak assumptions, we show that the flow of the Wasserstein distance between two solutions is contractive. Under more stringent conditions ensuring the regularity of the probabilistic solutions, we finally derive an explicit formula for the time derivative of the flow and we deduce the convergence of solutions to equilibrium.Comment: Stochastic partial differential equations: analysis and computations (2013) http://dx.doi.org/10.1007/s40072-013-0014-

Adhesion and volume constraints via nonlocal interactions determine cell organisation and migration profiles

The description of the cell spatial pattern and characteristic distances is fundamental in a wide range of physio-pathological biological phenomena, from morphogenesis to cancer growth. Discrete particle models are widely used in this field, since they are focused on the cell-level of abstraction and are able to preserve the identity of single individuals reproducing their behavior. In particular, a fundamental role in determining the usefulness and the realism of a particle mathematical approach is played by the choice of the intercellular pairwise interaction kernel and by the estimate of its parameters. The aim of the paper is to demonstrate how the concept of H-stability, deriving from statistical mechanics, can have important implications in this respect. For any given interaction kernel, it in fact allows to a priori predict the regions of the free parameter space that result in stable configurations of the system characterized by a finite and strictly positive minimal interparticle distance, which is fundamental when dealing with biological phenomena. The proposed analytical arguments are indeed able to restrict the range of possible variations of selected model coefficients, whose exact estimate however requires further investigations (e.g., fitting with empirical data), as illustrated in this paper by series of representative simulations dealing with cell colony reorganization, sorting phenomena and zebrafish embryonic development

Structure preserving schemes for mean-field equations of collective behavior

In this paper we consider the development of numerical schemes for mean-field equations describing the collective behavior of a large group of interacting agents. The schemes are based on a generalization of the classical Chang-Cooper approach and are capable to preserve the main structural properties of the systems, namely nonnegativity of the solution, physical conservation laws, entropy dissipation and stationary solutions. In particular, the methods here derived are second order accurate in transient regimes whereas they can reach arbitrary accuracy asymptotically for large times. Several examples are reported to show the generality of the approach.Comment: Proceedings of the XVI International Conference on Hyperbolic Problem

Single to double mill small noise transition via semi-Lagrangian finite volume methods

We show that double mills are more stable than single mills under stochastic perturbations in swarming dynamic models with basic attraction-repulsion mechanisms. In order to analyse this fact accurately, we will present a numerical technique for solving kinetic mean field equations for swarming dynamics. Numerical solutions of these equations for different sets of parameters will be presented and compared to microscopic and macroscopic results. As a consequence, we numerically observe a phase transition diagram in terms of the stochastic noise going from single to double mill for small stochasticity fading gradually to disordered states when the noise strength gets larger. This bifurcation diagram at the inhomogeneous kinetic level is shown by carefully computing the distribution function in velocity space

Explicit flock solutions for Quasi-Morse potentials

We consider interacting particle systems and their mean-field limits, which are frequently used to model collective aggregation and are known to demonstrate a rich variety of pattern formations. The interaction is based on a pairwise potential combining short-range repulsion and long-range attraction. We study particular solutions, that are referred to as flocks in the second-order models, for the specific choice of the Quasi-Morse interaction potential. Our main result is a rigorous analysis of continuous, compactly supported flock profiles for the biologically relevant parameter regime. Existence and uniqueness is proven for three space dimension, whilst existence is shown for the two-dimensional case. Furthermore, we numerically investigate additional Morse-like interactions to complete the understanding of this class of potentials.Comment: 26 page

Computation of power law equilibrium measures on balls of arbitrary dimension

We present a numerical approach for computing attractive-repulsive power law equilibrium measures in arbitrary dimension. We prove new recurrence relationships for radial Jacobi polynomials on d-dimensional ball domains, providing a substantial generalization of the work started in Gutleb et al. (Math Comput 9:2247–2281, 2022) for the one-dimensional case based on recurrence relationships of Riesz potentials on arbitrary dimensional balls. Among the attractive features of the numerical method are good efficiency due to recursively generated banded and approximately banded Riesz potential operators and computational complexity independent of the dimension d, in stark constrast to the widely used particle swarm simulation approaches for these problems which scale catastrophically with the dimension. We present several numerical experiments to showcase the accuracy and applicability of the method and discuss how our method compares with alternative numerical approaches and conjectured analytical solutions which exist for certain special cases. Finally, we discuss how our method can be used to explore the analytically poorly understood gap formation boundary to spherical shell support

ASYMPTOTIC FLOCKING DYNAMICS FOR THE KINETIC CUCKER-SMALE MODEL

In this paper, we analyse the asymptotic behavior of solutions of the continuous kinetic version of flocking by Cucker and Smale [16], which describes the collective behavior of an ensemble of organisms, animals or devices. This kinetic version introduced in [24] is here obtained starting from a Boltzmann-type equation. The large-time behavior of the distribution in phase space is subsequently studied by means of particle approximations and a stability property in distances between measures. A continuous analogue of the theorems of [16] is shown to hold for the solutions on the kinetic model. More precisely, the solutions will concentrate exponentially fast their velocity to their mean while in space they will converge towards a translational flocking solution

On the singularity formation and relaxation to equilibrium in 1D Fokker–Planck model with superlinear drift

We consider a class of Fokker–Planck equations with linear diffusion and superlinear drift enjoying a formal Wasserstein-like gradient flow structure with convex mobility function. In the drift-dominant regime, the equations have a finite critical mass above which the measure minimising the associated entropy functional displays a singular component. Our approach, which addresses the one-dimensional case, is based on a reformulation of the problem in terms of the pseudo-inverse distribution function. Motivated by the structure of the equation in the new variables, we establish a general framework for global-in-time existence, uniqueness and regularity of monotonic viscosity solutions to a class of nonlinear degenerate (resp. singular) parabolic equations, using as a key tool comparison principles and maximum arguments. We then focus on a specific equation and study in more detail the regularity and dynamics of solutions. In particular, blow-up behaviour, formation of condensates (i.e. Dirac measures at zero) and long-time asymptotics are investigated. As a consequence, in the mass-supercritical case, solutions will blow up in L∞ in finite time and—understood in a generalised, measure sense—they will eventually have condensate. We further show that the singular part of the measure solution does in general interact with the density and that condensates can be transient. The equations considered are motivated by a model for bosons introduced by Kaniadakis and Quarati (1994), which has a similar entropy structure and a critical mass if d ≥ 3